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We have created vortices in two-component Bose-Einstein 
condensates. The vortex state was created through a coherent 
process involving the spatial and temporal control of intercon- 
version between the two components. Using an interference 
technique, we map the phase of the vortex state to confirm 
that it possesses angular momentum. We can create vortices 
in either of the two components and have observed differences 
in the dynamics and stability. 



The concept of a vortex is at the center of our un- 
derstanding of superfluidity. A vortex is a topological 
feature of a superfluid - in a closed path around the 
vortex, the phase undergoes a 2tt winding and the su- 
perfluid flow is quantized. Following the experimental 
realization of a dilute atomic Bose-Einstein condensate 
(BEC) ||l|, much theoretical effort has been directed to- 
wards the formation and behavior of vortices in atomic 
BEC This paper presents the experimental realiza- 

tion and imaging of a vortex in BEC. We use the method 
proposed by Williams and Holland to create vortices 
in a two-component BEC. An interference technique is 
used to obtain phase images of the vortex state and con- 
firm the 2n phase winding required by the quantization 
condition. We have also carried out preliminary studies 
of the stability of the vortices. 

Vortices are created in superfluid helium by cooling a 
rotating bucket of helium through the superfluid transi- 
tion, and a vortex forms for each unit of angular momen- 
tum. This does not work for BEC because it is formed 
in a harmonic magnetic trap. When the condensate first 
forms it occupies a tiny cross-sectional area at the center 
of the trap and is too small to support a vortex. Even- 
tually, the condensate grows to sufficient size that it can 
support vortices, but the time scale for vortices to be 
generated in the vortex-free condensate due to coupling 
with the rotating environment is unknown, and may well 
be longer than the lifetime of the condensate. This is the 
potential difficulty with using an optical "stirring beam" 
or magnetic field distortion to rotate the cloud during 
condensation, as has been frequently proposed. Another 
proposal has been to use optical beams with appropriate 
topologies to "imprint" a phase on an existing conden- 
sate. This technique must drive the local density to zero 
at some point and then rely on uncertain dissipative pro- 
cesses for the condensate to relax into a vortex state. 



We have avoided these uncertainties by creating vor- 
tices using a coherent process that directly forms the de- 
sired vortex wave function via transitions between two 
internal spin states of ^^Rb. The two spin states, hence- 
forth |1) and |2), are separated by the ground-state hy- 
perfine splitting and can be simultaneously confined in 
identical and fully overlapping magnetic trap potentials. 
A two-photon microwave field induces transitions be- 
tween the states. As we have seen in previous experi- 
ments, this coupled two-component condensate is exempt 
from the topological rules governing single-component 
superfiuids Q - rules that make it difficult to implant 
a vortex within an existing condensate in a controlled 
manner. In the coupled system, we can directly create a 
|2) (or |1)) state wave function having a wide variety of 
shapes |^ out of a |1) (or |2)) ground-state wave func- 
tion by controlling the spatial and temporal dependence 
of the microwave-induced conversion of |1) into |2). 

We control the conversion by shifting the transition 
frequency using the AC Stark effect. A spatially inho- 
mogeneous and moveable optical field (a focused laser 
beam) provides the desired spatial and temporal control 
of the AC Stark shift. The vortex state is an axially sym- 
metric ring with a 27r phase winding around the vortex 
core where the local density is zero. To create a wave 
function with this spatial symmetry, the laser beam is 
rotated around the initial condensate as in Fig. pi. The 
desired spatial phase dependence is obtained by detuning 
the microwave frequency from the transition, and rotat- 
ing the laser beam at the appropriate frequency uj to 
make the coupling resonant. For large microwave detun- 
ings S, the necessary rotation frequency is simply S. For 
smaller detunings, the rotation frequency must be the 
effective Rabi frequency of the microwave transition [Q. 
As shown in Fig. 03 for large detunings, the energy reso- 
nance condition now means that atoms can only change 
internal state through the coupling of the time-varying 
perturbation, and are therefore obliged to obey any se- 
lection rule that the spatial symmetry of that perturba- 
tion might impose. The center of the condensate (at the 
axis of the beam rotation) feels no time-varying change, 
while regions near the circumference of the condensate 
feel a near-sinusoidal variation, with a phase delay equal 
to the azimuthal angle 9 around the circumference of the 
cloud. Williams and Holland show that this is precisely 
the geometry best suited to couple the condensate into a 
vortex state. It should be emphasized that it is not sim- 
ply the mechanical forces of the optical field that excite 
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the vortex: a laser beam rotating clockwise can produce 
clockwise or counterclockwise vortex circulation, depend- 
ing on the sign of the microwave detuning. 




FIG. 1. (a) A basic schematic of the technique used to 
create a vortex. An ofF-resonant laser provides a rotating 
gradient in the AC Stark shift across the condensate as a 
microwave drive of detuning S is applied, (b) A level diagram 
showing the microwave transition to very near the |2) state, 
and the modulation due to the laser rotation frequency that 
couples only to the angular momentum I — 1 state when 
uj = 5. In the figure, the energy splitting (< 1 Hz) between 
the I = 1 and I = states is exaggerated. 

In the absence of the microwave coupling field, the two 
states can be thought of as two distinguishable, interpen- 
etrating superfluids that interact with each other and 
with themselves via a mean-field repulsion proportional 
to the local densities. The interaction coefficients differ 
shghtly 11:^, so the |1) fluid has slight positive buoyancy 
with respect to the |2) fluid @. When the |1) fluid 
has a net angular momentum, it forms an equatorial ring 
around the central |2) fluid. The |1) fluid partially pene- 
trates the constant-phase |2) fluid, which creates a central 
potential barrier. Conversely, a |2) vortex forms a ring 
that tends to contract down into the |1) fluid. We use the 
overlap of the |1) and |2) fluids to image the phase pro- 
file of the vortex state via the interconversion interference 
technique that we introduced in |jll[] . In the presence of a 
near-resonant microwave field (and no perturbing optical 
field) , the two states interconvert at a rate sensitive to the 
local difference in the quantum phases of the two states. 
Thus the application of a resonant 7r/2 microwave pulse 
transforms the original two-fluid density distribution into 
a distribution that reflects the local phase difference, a 
"phase interferogram." Looking at the condensate both 
before and after the interconversion pulse provides im- 
ages of both the amplitude and phase of the vortex ring. 

The basic experimental setup for forming condensates 
and driving them between different spin states is the same 
as in Using laser cooling and trapping, followed by 
trapping in a TOP magnetic trap and evaporative cool- 
ing, we produce a condensate of typically ~ 8 x 10^ atoms 
in the |1) state {F — 1, mp — —1). We then adiabati- 



caly convert the trap to a spherically symmetric potential 
by reducing the quadrupole magnetic field gradient . 
This leaves us with a condensate 54 microns in diameter 
in a trap with oscillation frequencies of 7.8 ±0.1 Hz in 
the radial and axial directions for both spin states. In 
this trap, a |1) state condensate has a lifetime of 75 s 
and the |2) state {F — 2, mp ~ +1) about 1 s. Oscil- 
lating magnetic fields can then be pulsed on to drive the 
microwave transition between the |1) and |2) states. The 
power and oscillation frequency of the fields are adjusted 
to obtain the desired effective Rabi frequency for the |1) 
to 1 2) transition. To create a vortex in |2) we add a 10 
nW, 780 nm laser beam that has a waist of 180 fim and 
is detuned 0.8 GHz blue of resonant excitation of the |2) 
state. Using piezoelectric actuators we rotate the beam 
in a ~75 /im radius circle around the condensate at 100 
Hz. The following procedure is used to obtain the pre- 
cise location of the laser beam that is required to create 
a vortex. We first set the effective Rabi frequency to 100 
Hz by adjusting the detuning of the microwave frequency 
away from resonance in the presence of the light. Typi- 
cally the resonant Rabi frequency is 35 Hz and the total 
detuning about 94 Hz from the |1) to |2) transition. By 
making fine adjustments of the detuning we then opti- 
mize the amount of transfer (typically 50%) to the |2) 
state. We then adjust the center of rotation of the beam 
to obtain the most symmetric rings. After ^ 70 ms the 
vortex has been formed, and we turn off the laser beam 
and the microwave drive. 

We can take multiple images of a vortex both during 
and after formation using nondestructive state-selective 
phase contrast imaging |^,|l^. Rapid control of the mi- 
crowave power and frequency allows us to apply various 
pulse sequences to explore many different options for the 
creation, manipulation, and observation of a single vor- 
tex. For example, we can put the initial condensate into 
either the |1) or |2) state and then make a vortex in the 
|2) or |1) state respectively, and wc can obtain phase in- 
terferograms or quickly switch the internal state of the 
vortex at any time after the rotating laser beam is off. 
We can also watch the evolution of the vortex over time 
scales from milliseconds to seconds. All of these tech- 
niques are nondestructive to both the density and phase. 

In Fig. ||we show a detailed picture of the phase profile 
of a 1 2) vortex. To obtain this we first take a picture of 
the vortex (Fig. j^a) , then we apply a resonant microwave 
TT-pulse. Half-way through the 4 ms long 7r-pulse, we 
take a second image (Fig. ||b), and at the completion 
of the pulse we take a third image (Fig. ^) that shows 
the original density distribution of the interior |1) state. 
Normalizing by the density distributions of the vortex 
and interior states we obtain the phase image in Fig. ^ 
Jl^ . The figure dramatically shows the variation and 
continuity of the phase around the ring (Fig. |^) that 
are required by the quantization of angular momentum. 
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FIG. 2. Condensate images (100 /xm on a side) in which the imaging laser is tuned such that only |2) is visible, (a) The 
intial image of the vortex, (b) after n/2 of the interconversion pulse, and (c) after completion of the 7r-pulse. The vortex is now 
invisible (in |1)) and the interior fluid is imaged (in |2)). (d) The normalized difference in densities between the local average 
ring and interior densities ((a) and (c) respectively) and the phase interferogram density (b) for each corresponding point in 
the images. This is approximately the cosine of the local phase difference <j) between the vortex state and the interior state [14] . 
The values are shown only for those regions where the densities of each state were high enough to give adequate signal to noise 
for the phase reconstruction, (e) The radial average at each angle 6 around the ring is shown in (d). (The data is repeated 
after the azimuthal angle 2n to better show the continuity around the ring.) 
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In Fig. ^ are pictures of the time evolution of vortices 
to show their dynamics and stabihty. As expected, the 
dynamics of the |1) state vortices (Fig. |^ and Fig. ||b) 
are different from the |2) state vortices (Fig. ||c) due to 
the different scattering lengths. The |2) vortex ring sinks 
in toward the trap center (Fig. |^c), and then rebounds 
and apparently breaks up. This pattern is repeatedly ob- 
served in measurements of the |2) vortices. Conversely, 
the equilibrium position of a |1) state fluid is obtained 
by "floating" outside the |2) state fluid. The inner ra- 
dius of the |1) state ring shrinks slowly as the interior |2) 
fluid decays away with a ^ 1 s lifetime. A variety of ad- 
ditional behaviors has been seen for |1) vortices. Initial 
asymmetry is very sensitive to beam position and con- 
densate slosh. For small differences in the initial vortex 
state density distribution, asymmetric density distribu- 
tions sometimes develop and/or "heal" between 0.5 and 
1 s. 

The nonrotating |2) fluid within a |I) vortex is analo- 
gous to the defects that pin vortices in superconductors. 
This 1 2} "defect" can be removed quickly, by a properly 
tuned laser pulse, or allowed to decay slowly (as in Figs. 
^ and . Thus this system can be varied between two 
relevant physical limits. In the limit of a large repul- 
sive central potential (produced by a large amount of |2) 
fluid) the system most closely resembles quantized flow in 
a fixed, 3-D toroidal potential. The vortex core is pinned 
and its size is determined by the central potential. Figs. 
^ and Ijb are near this limit for t~ and 200 ms. In this 
case the density distribution of the vortex state is vul- 
nerable to instabilities due to the fact that the relative 
densities of |1) and |2) may evolve with relatively little 
energy cost so long as the total density remains constant 
pO| . In the opposite limit (a small amount of |2} fluid), 
the central potential is negligible and the size of the vor- 
tex core (in equilibrium) is determined entirely by its own 
centrifugal barrier. Figs. ^ and ^3 at t = 600 ms are 
evolving towards this limit. In Fig. we see that the 
interior |2) fluid is no longer pinning the vortex at 600 
ms, but rather is being dragged around by the precessing 
vortex. The |2) vortex provides an interesting mixture of 
the above limits. The |I) fluid "floats" to the outside so 
there is no pinning of the vortex core, but the tendency 
towards density instabilities remains. 

Expanded studies of stability issues are underway. We 
also expect to be able to observe interesting transitional 
behavior between these limiting cases. For example, it is 
a straightforward extension of our method to create an 
I = 2 vortex. In the presence of a strong pinning po- 
tential, I = 2 vortices should be stable, but in the weak 
potential limit, I = 2 vortices are predicted to sponta- 
neously bifurcate [Q. We gratefully acknowledge useful 
conversations with J. Williams and M. Holland. This 
work is supported by the ONR, NSF, and NIST. 
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FIG. 3. (a,b) Two separate instances of the free evolution 
of a |1) state vortex in the magnetic trap. It is stable over a 
time long compared to the trap oscillation period (128 ms). 
(c) The free evolution of a j2) vortex is much more dynamic. 
It is seen shrinking quickly into the invisible |1) fluid and 
rebounding into fragments. Each column is from a single run, 
where time t is referenced to the end of vortex creation {t is 
the same for each row). The |1) and |2) state images appear 
different due to different signs of the probe detuning. 
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